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Abstract. We present ESO-NTT low resolution spectroscopy of the neutron star X-ray transient Cen X-4 in quiescence over 
a complete orbital cycle. Our data reveal the presence of a K3-7 V companion which contributes 63% to the 5600-6900 A 
^ ' flux and orbits the neutron star with a velocity semi-amplitude of K2 = 145.8 ± 1.0 km s"'. This, combined with a previous 

fT^ ' determination of the inclination angle and mass ratio, yields a neutron star and companion mass of Mi = 1.5 ± 1.0 Mq and 

00 . M2 = 0.31 + 0.27 Mq, respectively. The mass donor is thus undermassive for the inferred spectral type indicating it is probably 

10 ' evolved, in agreement with previous studies. Doppler tomography of the Ho- line shows prominent emission located on the 

00 I companion and a slightly asymmetric accretion disc distribution similar to that seen in systems with precessing eccentric discs. 

Strong Ho- emission from the companion can be explained by X-ray irradiation from the primary. No evidence is found for a 
hot spot in Ha, whereas one is revealed via Doppler tomography of the Hel lines. This can be interpreted as the hot spot and 
outer regions of the disc being at a higher temperature than in other systems. 

Oh' Key words, accretion, accretion discs — binaries: close — star: individual (Cen X-4) — stars: neutron 
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c/3 1 . Introduction which the X-ray luminosity drops by some 5-6 orders of mag- 

^ ' nitude. Concerning the optical properties, SXTs in outburst are 

• In recent years coherent, fast pulsations have been discovered relatively bright and their emission is dominated by X-ray re- 

'O '''' ^"'"''^ Binaries (LMXBs) leading to the con- processing in the accretion disc surrounding the neutron star, 

elusion that they contain rapidly spinning (-2-5 ms), weakly quiescence their optical luminosity drops by as much as 6-7 

. magnetic neutron stars (see e.g. Chakrabarty 2004). X-ray light j^ag, thereby allowing the late type companion to be detected 
- - ' curves obtained with RXTE also revealed coherent oscillations ^^^^^^^ Weak, residual optical emission is still present in 

during Type I X-ray bursts (i.e. thermonuclear explosions on quiescence, usually at the ~ 10 - 30% level in V, that is usu- 

the neutron star surface) in fourteen LMXBs (e.g. van der Klis ^jjy ascribed to a dim, viscously - heated accretion disc (van 

2000). These findings directly confirmed evolutionary models Paradijs & McCHntock 1995) 
that link the LMXB neutron stars to those of millisecond ra- 
dio pulsars, the former being the progenitors of the latter (e.g. 

Bhattacharya & van den Heuvel 1991). The large luminosity variations of SXTs correspond to vast 

A large fraction of LMXBs accrete matter at very high rates changes in the accretion rate onto the collapsed objects. These 

(and therefore shine as bright X-ray sources) only sporadically, provide an opportunity to explore a range of accretion regimes 

Among these systems are Soft X-ray Transients (SXTs), host- onto the neutron star that are inaccessible to persistent sources, 

ing old neutron stars (for a review see Campana et al. 1998). While in outburst accretion onto the neutron star surface takes 

These systems alternate periods (weeks to months) of high X- place in a manner similar to that of persistent LMXBs, whereas 

ray luminosity with long (1-10 years) intervals of quiescence in when the accretion rate decreases, matter is halted at the mag- 

netospheric boundary by a centrifugal baiTier generated by the 

's^l offprint requests to: Paolo D'Avanzo, da- fast rotating neutron star magnetosphere (a.k. a. propeller effect; 

vanzo@merate.mi.astro.it Illarionov & Sunyaev 1975). If the mass inflow rate decreases 

* The results reported in this paper ai'e partially based on observa- further, the neutron star can enter the radio pulsar regime and 

tions carried out at ESO, La Silla, Chile (67.D-01 16). begin to expel the inflowing material from the companion star 
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(Shaham & Tavani 1991; Stella et al. 1994; Campana et al. 
1995;Burderietal. 2001). 

Some progress has recently been made in the study of the 
outburst decay and quiescent state of some of these systems 
(Asai et al. 1996), yet the interpretation of the quiescent opti- 
cal and X-ray observations is ambiguous. The quiescent X-ray 
spectrum is well fit by the sum of a soft component (usually a 
neutron star atmosphere model with temperature kT ~ 0.1-03 
keV) plus a hard tail power law with photon index in the range 
1-2. There are cases in which only one of the two components 
is detected. Different interpretations have been put forward. 
Menou et al. (1999) suggested the presence of an Advection 
Dominated Accretion Flow (ADAF) dominating the residual 
optical emission: some of the inflowing matter would leak 
through the magnetosphere and accrete onto the neutron star 
surface, giving rise to the soft X-ray component. Campana et al. 
(1998a, b, see also Campana & Stella 2000) suggested instead 
that the soft component arises from the cooling of the neutron 
star surface heated during the repeated outburst episodes (see 
also Brown et al. 1998; Rutledge et al. 1999) and the hard tail 
arises from shock emission powered by an active radio pulsar 
In this case (part of) the residual quiescent optical flux may 
come from the shock emission in the interaction between the 
radio pulsar wind and the matter inflowing from the compan- 
ion. 

We therefore decided to utilise Doppler tomography in 
order to study the kinematics and geometry of matter in 
such interacting binaries. This would allow us to derive two- 
dimensional velocity maps of the optical line-emitting regions, 
by combining spectra taken at different orbital phases (Marsh 
& Horne 1988). This method provided strong evidence for the 
presence of accretion discs in cataclysmic variables (CVs) and 
quiescent black hole binaries (e.g. Casares et al. 1997). In the 
case of the CV AE Aqr, Doppler tomograms revealed a highly 
unusual configuration, which supported the idea that the white 
dwarf expels most of the incoming material from the com- 
panion star and is thus in the propeller regime (Wynn et al. 
1997; Horne 1999). The potential of this technique has not yet 
been fully exploited for neutron star SXTs. Only recently has a 
Doppler map of the SXT Cen X-4 been obtained (Torres et al. 
2002), although this required observations acquired over a year 
In addition, Doppler maps of XTE J2123-058 in both outburst 
and quiescence have been produced by Hynes et al. (2001) and 
Casares et al. (2002), respectively. 

2. Cen X-4 

Bright X-ray outbursts were detected from Cen X-4 in 1969 
(Conner et al. 1969) and again in 1979 (Kaluzienski et al. 
1980), since when it has been quiescent. The companion star 
is an evolved ~ 0.7 Mq, K5-7 star (Shahbaz et al. 1993; Torres 
et al. 2002) making Cen X-4 one of the brightest SXTs in qui- 
escence with V=18.7 and (Ay - 0.3 mag). The stellar parame- 
ters have recently been improved by Gonzalez Hernandez et al. 
(2005) who found a surface temperature of T^ff = 4500 ± 100 
K and surface gravity of log^ = 3.9 ± 0.3. A 15.1 hr orbital 
period was determined from the sinusoidal variation of the op- 
tical light curve (Cowley et al. 1988; Chevalier et al. 1989; 



Table 1. Radial velocity parameters of Cen X-4. 



Template Star 


Spectral Type 


y 




K2 




A. red 






(km s" 




(km s" 






HR 5536 


K3 III 


203.6 + 


2.1 


147.7 ± 


1.0 


1.8 


HD 184467 


K2 V 


191.2 ± 


2.1 


147.6 ± 


0.9 


1.9 


HD 29697 


K3 V 


209 ± 


10 


146.1 + 


1.0 


1.2 


HD 154712A 


K4 V 


194.2 ± 


2.1 


146.2 + 


0.8 


1.7 


61 Cyg A 


K5 V 


196.4 + 


1.0 


146.4 + 


0.8 


1.9 


61 Cyg B 


K7 V 


196.0 ± 


1.1 


145.8 ± 


1.0 


1.1 



McClintock & RemiUard 1990). These characteristics make 
Cen X-4 an exceUent candidate for Doppler tomography stud- 
ies in quiescence. 

The optical spectrum shows clear residual disk emission 
after subtraction of the companion's spectrum. Assuming a 
smooth spectral energy distribution, it contributes 80% in B, 
30% in V, 25% in R and 10% in 1 (Shahbaz, Nay lor & Charles 
1993). Strong Ha {EW ~ 40 A), Uj3 (EW ~ 20 A) and Hy 
(EW ~ 10 A) emission is also present. Although variable, these 
lines are always present in the optical spectra taken at different 
epochs. For instance, van Paradijs et al. (1989) reported equiv- 
alent widths of 43, 33 and 28 A, respectively. 

HST observations revealed a flat UV spectrum (in a vFy 
vs. V representation McClintock & RemiUard 2000), which 
connected smoothly with the hard tail X-ray power law (e.g. 
Campana & Stella 2000), suggestive of a shock spectrum. This 
is at variance with observations of the black hole candidate 
A0620-00. 



3. Observational data and reduction 

Cen X-4 spectra were taken on 26-28 May 2001 with the ESO 
3.5m New Technology Telescope (NTT) at La Silla. We used 
the ESO Multi Mode Instrument (EMMI) with a Tektronix 
CCD of 2048 X 2048 pixels yielding a field of view of 9. 1 ' x 8 .6' 
at a resolution of 0.27"/pixel. All nights were photometric with 
seeing in the 0.5 - 1" range. We used a 1" slit and grating 7 to 
cover the Ha region, centered at 6250 A, and with a dispersion 
of0.66A/pix. 

The data were reduced using standard ESO-MIDAS pro- 
cedures for bias subtraction, flat-field correction and cosmic 
ray removal. AU spectra were sky-subtracted and corrected for 
atmospheric extinction. Helium-argon lamp images were ob- 
tained for wavelength caUbration during daytime and with the 
telescope verticaUy parked. The wavelength scale was then de- 
rived through third-order polynomial fits to 26 lines, result- 
ing in an rms scatter of < 0.06 A. Instrumental flexure during 
our observations were then accounted for using atmospheric 
emission lines in the sky spectra. Spectral resolution was about 
~ 2.2 pixels, or ~ 70 km s~^ in velocity. 

Observations were carried out over three consecutive nights 
so as to completely cover the 15.1 hr orbital cycle of Cen X- 

4. We obtained 30 spectra, each 30 min long. This is the first 
time that an extensive and contiguous orbital coverage has been 
obtained for Cen X-4. 
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Reference J P 

(kms-') (kms-') (d) 

Cowley et al. 1988 234 ± 8 152 ± 10 '- 

McClintock & Remillard 1990 137 + 17 146 ± 12 0.629063 

Torres et al. 2002 184.0 ± 1.9 150.0 ± 1.1 0.6290496 + 0.0000021 

This paper 196.0 ±1.1 145.8 ± 1.0 0.630 ± 0.001 

The velocity values reported here have been obtained by fixing P to the value given by Torres et al. (2002). We also repeated the fit leaving all 
the parameters free, obtaining the value given here. 



4. Cen X-4 properties 

4. 1. Radial velocity curve 

Radial velocities were measured from our 30 spectra through 
cross-correlation with late-type template stars. Prior to cross- 
correlating, the spectra were rebinned to a uniform velocity 
scale. Cross-correlation was performed in the range 5930-6750 
A, after masking out strong night-sky lines due to OI 6300 A, 
the main atmospheric and interstellar features, emission lines 
from the accretion disc {}\a, Hel /15876 and /J6678) and our 
spectral extremes, resulting in an effective wavelength range 
of 6380-6520 A. We then performed a (least-squares) sine fit 
to our velocity data. We adopted the parameters obtained with 
the K7 V template (see Table as they provided the best fit 
in terms of reduced giving K2 - 145.8 + 1.0 km s"', and 
y = 196.0 ± 1.1 km s"', which are, respectively, the radial ve- 
locity amplitude and systemic velocity of Cen X-4. 

The dynamical parameters obtained here are given in Table 
IS together with the results of Cowley et al. (1988), McClintock 
& Remillard (1990) and Torres et al. (2002). Our value for K2 
is slightly smaller than that of ToiTes et al. (2002) and is statis- 
tically inconsistent. However, we note that the data of Torres 
et al. are spread over several years, whereas ours were col- 
lected over three consecutive nights. Our result is instead con- 
sistent with those of Cowley et al. (1988) and McClintock & 
Remillard (1990). Curiously, a much larger variation is found 
in the systemic velocity y, as has also been noted by ToiTes et 
al. (2002). We obtain a y intermediate to the values reported 
earUer. We do not have an explanation for this discrepancy, 
which in principle could be due to the use of an erroneous he- 
liocentric radial velocity for the template stars. Our value for P 
has a correspondingly large uncertainty as a result of our rela- 
tively short observing baseline with respect to, e.g. Torres et al. 
(2002). 

The radial velocity curve, after folding with our best 
ephemeris, is shown in Fig.[n It is highly sinusoidal (x^^, = 1.1) 
despite the point at phase 0.75. This point is not due to instru- 
mental problems (such as telescope flexure) and we also note 
that a similar effect is seen in the radial velocity curve of Cen 
X-4 measured by Torres et al. (2002); however there is no clear 
explanation for such behaviour. 

4.2. Spectral classification 

We have produced an average spectrum of Cen X-4 in the rest 
frame of the companion star after Doppler shifting every indi- 



vidual spectrum using our orbital solution and the orbital pe- 
riod of Torres et al. (2002). The spectrum shows prominent 
Ha, Hel /15876, /16678 emission as well as absorption fea- 
tures typical of late type stars (e.g. TiO /16161, Fell /16361 
and Fell /16494, Fig. |3- A spectral classification can be de- 
rived by taking different templates degraded appropriately with 
a Gaussian bandpass to match the broadening of our NTT ob- 
servations. These templates were then multiplied by a factor 
< / < 1, representing the fractional contribution of light 
from the secondary star and subtracted from the target average 
(after masking interstellar, night-sky and emission lines and re- 
binning to a uniform velocity scale). The best fit is obtained 
for templates K3 V-K5 V and / = 0.59 - 0.62, although a K7 
III also provides an acceptable fit (see Table |3}. More impor- 
tantly, in the Doppler-corrected average, Li I /16708 absorption 
becomes prominent. The detection of a strong Li /16708 ab- 
sorption feature in such a late-type companion would normally 
be unexpected since the star's initial Li I content should be 
rapidly depleted by a contamination of convective mixing and 
mass transfer to the compact object. Such unusual abundances 
of Li I have already been found by other authors (e.g. Martin 
et al. 1994, Torres et al. 2002). This significant and anoma- 
lous abundance of Li I in Cen X-4 (and other transients) finds 
various explanations in the literature, such as synthesis in a su- 
pernova explosion of the compact primary's progenitor, or a-a 
reactions during the repeated strong outbursts that characterize 
transient X-ray binaries (Martin et al. 1994), a relatively recent 
formation of the system (Gonzalez Hernandez et al. 2005), or 
an effect related to the tidally-locked rotation of the two stars, 
which leads to a slower lithium destruction rate (Maccarone et 
al. 2005). 

An alternative approach to the spectral classification con- 
sists of comparing the ratio of the equivalent width of the ab- 
sorption lines in Cen X-4 with different templates. We consider 
the ratios ^6361(FeII)M6161 (TiO) and ^6494(FeII)//16161 
(TiO), which also classify the secondary of Cen X-4 as a K5 
V-K7 V star (see Tablel^. We therefore conclude that the most 
likely spectral type for Cen X-4 is a K3-7 V star. 

4.3. Orbital variation of the companion star 
contribution 

Taking advantage of the good signal to noise of each individual 
spectrum we tried to estimate the factor / for all the spectra, 
allowing us to monitor its variability with orbital phase. It is 
interesting to note that / is not constant but is modulated at 
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Fig. 1. Radial velocity curve for Cen X-4 computed using a K7 V template. Observations were performed on 26 (triangles), 27 
(dots) and 28 (stars) May 2001. Two phases are shown for clarity. The best sine-wave fit is also shown (see Section 4.1). 



Table 3. Cen X-4 spectral classification by direct fitting 



Template Spectral Type / x v 



(117d.o.f.) 



HR5536 


K3 III 


0.49 ± 0.01 


2.29 


HR6159 


K4III 


0.44 ± 0.01 


2.40 


HR 5622 


K5 III 


0.53 ±0.01 


2.30 


HD 184467 


K7 III 


0.84 + 0.05 


1.30 


HD 29697 


K3 V 


0.62 ± 0.04 


1.29 


HD 1547 12A 


K4 V 


0.65 ± 0.04 


1.32 


61 Cyg A 


K5 V 


0.63 ± 0.04 


1.30 


61 CygB 


K7 V 


0.59 ± 0.04 


1.39 


Table 4. Line ratios for Cen X-4 and template stars. 


Template 


Spectral Type 


6361 (Fell)/ 


6494(FeII)/ 






6161 (TiO) 


6161 (TiO) 


Cen X-4 




0.41 ± 0.06 


0.52 ± 0.06 


HR 5536 


K3 III 


1.43 


1.52 


HR6159 


K4III 


1.59 


1.88 


HR 5622 


K5 III 


0.76 


0.86 


HD 184467 


K7III 


0.68 


0.78 


HD 29697 


K3 V 


0.18 


0.44 


HD 1547 12A 


K4 V 


0.12 


0.51 


61 Cyg A 


K5 V 


0.33 


0.52 


61 CygB 


K7 V 


0.50 


0.56 



the orbital period (see Fig.|3}. Variations of / with phase show 
two minima at phase (i.e. at superior conjunction) and 0.5 



and two nearly equal maxima at phase 0.25 and 0.75. The lat- 
ter occur when the observer sees the elongated star with max- 
imum projected surface area. This behaviour follows the clas- 
sical ellipsoidal modulation previously observed in optical and 
IR lightcurves (e.g. McClintock & Remillard 1990; Shahbaz 
et al. 1993), and strongly suggests that it is dominated by the 
visibility of the tidally distorted companion. However, we note 
that gravity darkening effects, expected from the classical ellip- 
soidal modulation, are not sufficient to explain the large depth 
difference between phase and 0.5 observed in our lightcurve. 
Such a large difference may require the presence of an addi- 
tional effect, such as the heating of the inner hemisphere of the 
companion by X-ray irradiation from the compact object. We 
give further evidence of such a phenomenon in Section 5, in the 
light of our Doppler tomography results. 

4.4. System parameters 

If we assume that the secondary star is syncronized with the 
binary motion and fills its Roche lobe, we can determine the 
mass ratio (q - M2/M1) through the expression (see e.g. Wade 
& Home 1988) 

V sin / = 0A62K2q^^^(l + qf'\ 

Taking the secondary star's rotational broadening as mea- 
sured by Torres et al. (2002) (v sin/ = 43+6 km s"'), and 
combined with our value of K2, we derive q - 0.18 + 0.06. 
Combining our value of K2 with the precise orbital period re- 
ported by Torres et al. (2002), we derive a mass function of 



D'Avanzo et al.: Doppler tomography of Cen X-4 



5 




Fig. 2. Doppler-corrected average of the 30 spectra of Cen X-4 in quiescence that were obtained and analyzed in the present study 
(top spectrum) compared with a K5 V template (61 Cyg A). Arbitrary vertical offsets have been applied for clarity. 



f{M) - 0.201 + 0.004 Mq. In addition, if we take again the 
V sin/ value measured by Torres et al. (2002) and the con- 
straint on / from the IR observations of Shahbaz et al. (1993, 
i.e. 31° < / < 54° at 90% confidence level), we can constrain 
the masses of the two stars: 0.47 < M^s < 2.31 Mq for the 
neutron star and 0.06 < M2 < 0.55 for the companion. M2 is 
clearly undermassive for the inferred spectral type as discussed 
by Shahbaz et al. 



5. Doppler tomography 

The three emission lines seen in the Cen X-4 spectra are 
double-peaked, testifying to the presence of an axisymmet- 
ric emission region around the compact object. The line pro- 
files vary with orbital phase and there is evidence for different 
emitting components. Algorithms have been developed to re- 
construct the velocity field of the line emitting region, i.e. the 
Doppler tomography technique (Marsh & Horne 1988). The 
continua were normalized and subtracted so as to give a pure 
emission line profile which was subsequently rebinned onto 
a constant velocity scale of 30.5 km s ' pixel The result- 
ing Doppler maps are presented in Fig. |3 in velocity coor- 
dinates, together with the trailed spectra. The Roche lobe of 
the companion and the theoretical path of the gas stream for a 
K2 - 145.8 km s"' and q - 0.18 system are superposed over 
each map, in steps of 0. 1 R/,, . 



5.1. Ha tomography 

The Ha emission line has a mean equivalent width (EW) of 
45.3 + 1.2 A. The Ha map is dominated by a remarkably in- 
tense spot, coincident in phase and velocity with the compan- 
ion star. This may indicate either chromospheric activity in the 
synchronized (and hence rapidly rotating) late-type companion 
or irradiation from the compact object. There is also evidence 
for a ring-like structure, a signature of a classic accretion disc. 
We note that the brightness distribution of the disc is slightly 
asymmetric, with enhanced emission at the negative V,, quad- 
rants. This is somewhat reminiscent of the emissivity distri- 
bution observed in SW Sex stars, a group of persistent cata- 
clysmic variables where exotic scenarios (namely magnetic ac- 
cretion, disc overflow and propeller mechanisms) have been in- 
voked to explain their peculiar emission line profiles. The same 
phenomenology is sometimes observed in LMXBs both dur- 
ing outburst and quiescence (e.g. XTE J2123-058: Hynes et 
al. 2001; Casares et al. 2002). We note, however, that one of 
the signatures of SW Sex stars is the presence of low-velocity 
emission in their Doppler maps (e.g. Wynn et al. 1997; Welsh, 
Horne & Gomer 1998). Although our Ha tomogram shows an 
asymmetric distribution, there is no evidence for emitting struc- 
tures at velocities lower than ~200 km s"' and hence a SW Sex 
scenario does not seems a likely model. Alternatively, asym- 
metric disc emissivity can be produced by spiral waves (trig- 
gered by tidal interaction with the companion star, see Steeghs, 
Harlaftis &. Horne 1998) or superhumps caused by eccentric 
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Fig. 3. Variation of fractional contribution / of the companion star with orbital phase. Observations were performed on 26 
(triangles), 27 (dots) and 28 (stars) May 2001. Two cycles are shown for clarity. 



precessing discs (Foulkes et al. 2004). Finally, direct irradia- 
tion of the outer disc rim can also excite enhanced emission. In 
this latter case, our asymmetric Ha tomogram may indicate a 
vertical irradiated structure around orbital phases 0.3-0.5, pos- 
sibly connected with a shock front. 

On the other hand, we see no sign of emission from the 
expected hot spot location. This pattern has also been observed 
in other quiescent (black hole) X-ray transients, i.e. Nova Mus. 
91 (Casares et al. 1997), Nova Oph. 77 (Harlaftis et al 1997), 
XTE J1118H-480 (Ton-es et al. 2005). 

In order to further investigate the possible origins of the Ha 
emission coming from the secondary, we extracted the narrow 
Ha component from our Doppler map. To do this we subtracted 
from our map a simulated axisymmetric image centered on the 
compact object (Vj - 0, V,. - -K2/q), obtaining a map where 
the emission arising from the secondary becomes clearer (Fig. 
refdopsub, first panel). Then, projecting this new Doppler im- 
age into the observed orbital phases and adding all the residu- 
als together we obtained a spectrum which represents the Ha 
emission from the secondary. Following Casares et al. (1997) 
and Torres et al. (2002) we estimated an equivalent width for 
the narrow line component of Ha of 4.4 + 0.5 A, taking into 
account the veiling from the accretion disc. This is quite a high 
value of Ha emission for such systems. Using the relationship 
given by Soderblom et al. (1993): 



logFH„ = logEW(H^) + 0.n3{B-V)l-l.lSS{B-V)o+7AS7 



we convert this value into an Ha line surface flux at the star 
(FHa)- Using (B - V)q = 1.25 mag (Chevalier et al. 1989) we 
obtain Fho = 7 x 10^ erg cm"^ s"'. 

A chromospherically active secondary is a possible expla- 
nation for this Ha emission. We can compare our estimates 
with the results reported by Soderblom et al. (1993) for rapidly 
rotating dwarfs in the Pleiades. From Table 6 of that paper we 
find that the Ha emission of the Cen X-4 companion is higher 
in strength than that observed for all the K dwarfs in their sam- 
ple. Moreover, we note that Torres et al. (2002) found for the 
narrow Une component of Ha an EW of ~ 2.3and a veiling fac- 
tor (1 -/) of 0.25 (where / is the fractional contribution of light 
from the secondary star) while our higher value of EW corre- 
sponds to a higher value of (1 - /) (0.37, see TableO, showing 
that the veiling and Ha EW are correlated. This is not expected 
in a chromospheric Ha emission scenario. 

X-ray irradiation by the compact object is an alternative ex- 
planation for the Ha emission from the secondary. The intrin- 
sic X-ray luminosity in the 0.5-10 keV band was estimated by 
Campana et al. (2004) to be L ^ 4 x 10^^ erg s"' . The X-ray ir- 
radiation at the secondary star is then Fx - Lxli^na^) - 5x10^ 
erg cm"^ s"', for an orbital separation a - 3.6 Rq, that can 
power the Ha emission if at least 1% of the incident X-ray flux 
is reprocessed to Ha photons. 

Following Hynes et al. (2002), a rough estimate of the frac- 
tion of the X-ray luminosity converted to Ha photons can be 
obtained assuming that the observed Ha emission is due com- 
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Fig. 4. Doppler maps of the Ha, Hel /15876 and Hel /16678 lines. The trailed spectra are shown in the upper panels and, for 
clarity, the same cycle has been plotted twice. The Roche lobe of the secondary and the predicted velocity of the gas stream are 
plotted for K2 - 145.8 km s ' and q = 0.18. Distances are marked along the curve in open circles and are in multiples of Q.lRn. 
The center of mass of the system and the neutron star position are denoted by crosses. 



pletely to X-ray reprocessing by the companion: 



Lho - fifiLx 



where /i is the fraction of the high-energy emission Lx which is 
intercepted by the companion (i.e. the solid angle of the com- 
panion) and /2 is the fraction of input energy emitted in Ha 
that we assume equal to 0.3 (Osterbrock 1987). The result is 
Lfja - 5 X lO^^Lx- This is in agreement (within a factor of 2) 
with our estimate of the fraction Lx converted to Ha photons 
and means that the incident X-ray flux would be sufficient to 
power the observed Ha emission. The correlation between the 
Ha EW and the veiling factor (1 - /) finds a natural explanation 
in this scenario. 



5.2. He line tomography 

Hel /15876 and Hel /16678 Doppler maps are the first to be ob- 
tained for a neutron star system. Hel maps, as for Ha, show a 
clear accretion disc pattern, with a hint of a hot-spot compo- 
nent, particularly intense in Hel /J6678. It is also apparent from 
the Doppler images that the disc extends to higher velocities 
in the Hel lines than in Ha, with approximate radii at ~300 
km s"' and ~200 km s"' respectively. This indicates that Hel 
emission arises from the inner regions of the disc, whereas Ha 
emission is more concentrated in the outer disc regions. This 
is as expected for viscously heated discs, where Te/f decreases 
outwards. 

Also for these lines, we subtracted the disc contribution 
to each map by simulating an axisymmetric image (Fig. |5}- 
After the subtraction, both Hel maps show emission in a 
region consistent with the theoretical position of the hot-spot. 
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-500 500 
Velocity (km/s) 

Fig. 5. Doppler maps of the Ha, Hel /15876 and Hel /16678 
lines (from top to bottom) with the axisymmetric part of the 
emission component subtracted. Plot symbols are the same as 
in figlH 



but only the Hel A5816 map shows emission coming from 
the secondary. These results put in a different perspective 
the earlier suggestions of a parallel with SW Sex systems. 
In particular, since the line energy increases from Ha to 
Hel (from 13.6 eV to 24.6 eV), it suggests that the hot spot 
temperature is high enough to emit fewer Ha photons but more 
He photons. Another indication of the different temperature 
of the emitting region comes from the measure of the position 
along the Vy axis of the emission arising from the companion 
in the Ha and Hel /15876 maps after the subtraction of the disc 
contribution: 



Ha: (Vy) = (97.7+ 11.8) km s-i 
HeI/15876: (Vy) = (122.8 + 11.8) km s^^ 

In a simple model where the accretion disc absorbs lower en- 
ergy photons more efficiently, we would expect lower excita- 
tion lines to be formed higher up in the Roche lobe and show 
higher velocities. A similar behaviour has been found in IP Peg 
by Harlaftis (1999). This scenario, however, cannot explain our 
different observed velocities. In fact, as we can see from our 
measurements, there is an indication that the centroid of Hel 
/i5876 is more shifted along the Vy direction than that of Ha. In 
any case, we can say that there is a region around ~ 80 km s"' 
from Li that may mainly be clear of emission due to the shadow 
of the disc, while the portion of the surface of the companion 
star which is directly irradiated by the neutron star shows clear 
Ha and Hel emission. 

This might indicate that the Ha emission from the sec- 
ondary comes mainly from a region near the Li point, where 
the radiation coming from the compact object is mostly oc- 
culted by the disc. Accordingly, the Hel /i5876 emission might 
come from a higher velocity region such as the portion of the 
surface of the companion star which is directly irradiated (and 
hence hotter) by the neutron star (i.e. it is not occulted by the 
disc). 

6. Conclusions 

Cen X-4 is the brightest (and nearest) quiescent SXT in the op- 
tical, leading to its being extensively studied in recent years. It 
is one of only a handful of SXTs for which Ha Doppler tomog- 
raphy has been obtained (in quiescence). Here we presented 
some refinements of the companion's characteristics, classify- 
ing it as an (undermassive) K3-7 V and improving on the sys- 
tem parameters. We produced the first light curve of the com- 
panion star contribution as a function of orbital phase, confirm- 
ing the shape and the model of optical and IR light curves of 
Cen X-4 (McCHntock & Remillard 1990; Shahbaz et al. 1993) 
and that the companion (almost) fills its Roche lobe. 

Doppler tomography represents an ideal tool to test and 
study the kinematics and geometry of the emission line regions 
of quiescent SXTs. Our Ha Doppler map shows significant 
emission at the Li point but no emission at the predicted lo- 
cation of a hot spot. The entire observed Ha emission arising 
from the companion seems to be powered by the quiescent X- 
ray luminosity of Cen X-4 (Campana et al. 2004), suggesting 
that the companion is irradiated by the compact object. 

The Ha Doppler map is not perfectly symmetric but shows 
enhanced emission at negative y velocities (see Fig.0}, possibly 
indicating superhumps caused by eccentric precessing discs. A 
similar pattern is also visible in the Ha Doppler map of Cen 
X-4 obtained by Torres et al. (2002). Hel A5816 and A6618 
maps add further detail to this picture. In particular, a hot spot 
component appears with a centroid at Hel /15876 that is more 
shifted along the Vy direction than that of Ha. These are in- 
dications for a temperature variation, with a hot spot able to 
emit higher energy lines and a cooler star emitting Ha and Hel 
/15876 but not Hel /16678. This can be obtained through irra- 
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diation of the disc and hot spot region and shielding of the Li 
point. 
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